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We propose a realistic scheme to create motional entangled states of a few bosonic atoms. It can 
experimentally be realized with a gas of ultra cold bosonic atoms trapped in a deep optical lattice 
potential. By simultaneously deforming and rotating the trapping potential on each lattice site it is 
feasible to adiabatically create a variety of entangled states on each lattice well. We fully address the 
case of Ai" = 2 and A*" = 4 atoms per well and identify a sequence of fractional quantum Hall states: 
the PfafBan state, the 1/2-Laughlin quasiparticle and the 1/2-Laughlin state. Exact knowledge of 
the spectrum has allowed us to design adiabatic paths to these states, with all times and parameters 
well within the reach of current experimental setups. We further discuss the detection of these 
states by measuring different properties as their density profile, angular momentum or correlation 
functions. 
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I. INTRODUCTION 

The creation of highly entangled multiparticle states is 
one of the most challenging goals of modern experimental 
quantum mechanics. In this respect atomic systems of- 
fer a very promising arena in which entangled states can 
be created and manipulated with a high degree of con- 
trol. The experimental difficulty increases, however, with 
the number of particles that are to be entangled, since 
the system becomes then more sensitive to decoherence. 
Starting with a small number of particles as a first step, 
important achievements have been already obtained in 
the creation of atomic entangled states. For example, in 
recent experiments with trapped ions, entangled states 
of up to four ions have been demonstrated 0. Moreover, 
in experiments with neutral bosonic atoms in optical lat- 
tices Bell-type states have been created by accurately 
controlling the interactions between neighbouring atoms 
As a typical feature of most of the experimentally 
realized entangled states, atoms get entangled through 
their internal degrees of freedom, keeping separable their 
motional part. 

In this article we develop a scheme to create motional 
entangled states of a small number of atoms in an actual 
experimental setup with an optical lattice These 
states are a sequence of fractional quantum Hall (FQH) 
states, analogous to the ones that appear in the context of 
fractional quantum Hall effect 0. In contrast to typical 
atomic entangled states, the particles are here entangled 
in real space, and not in internal space. This peculiarity 
makes them specially interesting, for it represents a novel 
nature of atomic entanglement. 

The possibility of creating FQH atomic states as the 
Laughlin state by rapidly rotating the trap confining the 
atoms has been discussed in several theoretical works 
However, experiments dealing with typically large 
number of particles have not yet succeeded in reaching 
these states. Here, we fully address the case of a small 
number of particles and design a realistic way of entan- 
gling them into FQH states. The experimental setup 



that we have in mind corresponds to a situation in which 
a Bose-Einstein condensate is loaded in a deep optical 
lattice. When the lattice depth is very large tunneling 
between different sites is strongly suppressed and the sys- 
tem can be treated as a lattice of independent wells, each 
of them with a small number of particles. By indepen- 
dently rotating each of these 3D wells Q the lowest Lan- 
dau level (LLL) regime can be achieved for each copy. 
We have studied the problem exactly within the LLL for 
N — 2 and iV = 4 particles per well. We have identified a 
sequence of highly entangled stable ground states, which 
are the Pfaffian state jlflj, the 1/2-Laughlin quasiparticle 
|Tll | and the 1/2-Laughhn state [ill- The 1/2-Laughhn 
quasiparticle state (which had never been identified be- 
fore in an atomic system) is particularly interesting. It 
is the counterpart of the 1/2-Laughlin quasihole found 
in and contains a 1/2-anyon. Driving the system into 
these strongly correlated states is, however, not trivial. 
By simply increasing the frequency of rotation the sys- 
tem will stay in a trivial non-entangled state with angu- 
lar momentum zero. Exact knowledge of the spectrum 
of the system has allowed us to design adiabiatic paths 
to these states by simultaneously rotating and deforming 
each of the wells. All parameteres and evolution times 
lie well within the reach of present experimental setups. 
We further discuss how to detect these entangled states 
by measuring different properties as their density profile, 
angular momentum or correlation functions. In particu- 
lar, we propose a novel technique to measure the density- 
density correlation function of these strongly correlated 
states. Even though the number of atoms per well is 
small, the lattice setup allows to have multiple copies 
of the system, so that the experimental signal is highly 
enhanced. 

We point out that our findings also show that adia- 
batically achieving FQH states for rapidly rotating traps 
with a large number of particles turns out to be very 
challenging, since the relevant experimental parameters 
scale linearly with the number of particles. Nevertheless, 
we hope that our results can shed some light on the prob- 
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lems that these current experiments are dealing with, and 
even may pave the way to new methods of achieving FQH 
multiparticle entangled states. 



II. IDENTIFICATION OF ENTANGLED 
STATES 

We consider a system of bosonic atoms loaded in a 3D 
optical lattice. We assume a commensurate filling of N 
atoms per lattice site [l^ , and a large value of the lattice 
depth Vq/Ej^^ > 1, where = fi^k'^/2M is the recoil 
energy, k is the wave vector of the laser lattice light, 
and M the atomic mass. In this limit the lattice can by 
treated as a system of independent 3D harmonic wells, 
each of them having N atoms and a trapping frequency 
w « VVoEr. 

Let us rotate each of these 3D harmonic wells around 
the direction X3 with frequency f2. We will identify a 
sequence of motional entangled ground states of the N 
atoms that appear as the frequency f2 is increased. We 
will assume the limit of rapid rotation In this case 
the motion in the ^3 direction is frozen, and the motion 
in the plane of rotation xi,X2 is restricted to the LLL. 
Note that in order to project the system onto the LLL 
we do not need to start with a 2D configuration (as it is 
the case in previous proposals since the fast rotation 
itself restricts the motion in the direction of the rotation 
to zero point oscillations. The system is then governed by 
a two dimensional effective Hamiltonian, which written 
in units of hiu has the form: 

H = {l-n/Lu)L + 2T: rjV, (1) 

where L = J2m=o''^'^m'^rn is the angular mo- 
mentum operator in the X3 direction, and V — 
Emi,m2,m3,m4^mi%"2'"mi"m2«'"3am4 IS the interaction 
operator. Here the bosonic operator aj„(a„i) create (ani- 
hilate) an atom in the state \m) of the LLL with well 
defined ^3 component of the angular momentum m. The 
wave functions of the LLL in complex coordinates read 

^^{z) = {z\m) = -jl=- z™e-l^l^/2 , (2) 

where z = {xi+ix2) / I = \/ h/Muj^ and m = 0, 1, . . . 00. 
Assuming contact interactions between the atoms the in- 
teraction coefficients are: 

mi,m2 2™i+™V"ii!m2!m3!m4! ' ^ ' 

In Hamiltonian |^ we have introduced the important in- 
teraction parameter 77 = ^J^pKasj with Os the 3D scat- 
tering length. Analytical calculations for scattering po- 
tentials of finite size ap have confirmed that the pseudo- 
potential approximation is also valid for tight traps with 
Os <C € as long as oq ^ £ is fulfilled [l^ . 



A. N=2 

First we consider the case of two particles per lattice 
well, which can be solved analytically. The Hamiltonian 
is diagonal in the states Im^, Wcm) of well defined rel- 
ative (jTir) and center of mass (mem) angular momentum: 

a = ^ £^m,,m,„|"lr,nT-cm)(?7lr,m-cm|, (4) 
with Era^,ra^^ = Sm^fiT] + (1 - ^/Lu){mr + mem)- We 

note that due to the restriction to s-wave scattering, only 
particles with zero relative angular momentum feel the 
interaction energy. It follows that for fi/w < 1 — 77/2 
the ground state of the system is |0, 0) (with total angu- 
lar momentum L = 0), which is not entangled, whereas 
for n/uj > 1-77/2 the state |2,0) (with L = 2) be- 
comes energetically favourable. This state, {zi,Z2\2,Q) cx 
{zi — 22)^ e~'^^' /^e"'^^' is clearly entangled since it 
cannot be written as a product of two single particle wave 
functions. It is the Laughlin state \iPl) for two particles 
at filling factor 1^ — 1/2 |ll|. In order to quantify the en- 
tanglement of this state we write it in the basis of states 
|mi77i2) with well defined single-particle angular momen- 
tum. Then the Laughlin state takes the form of a pure 
two qutrit state: [7/;^) = i (|02) + |20)) - ^|11)- This 
is already the Schmidt decoinposition of the state, and 
the entropy of entanglement [l6j can immediately calcu- 
lated to be E{\iPl)) = 1-5. This value is close to log2 3, 
corresponding to a maximally entangled pure two qutrit 
state. 



B. N=3 

The case of three particles per lattice well is very sim- 
iliar to the situation for N = 2. The 1/2-Laughlin state 
(L = 6) emerges as ground state after an intermediate 
state with odd angular momentum L = 3. As we will ex- 
plain in the next section, ground states with odd angular 
momentum cannot be reached using our proposal. Hence 
we now focus on a setup with four particles per lattice 
well, for which an interesting sequence of prominent FQH 
states arises. 



C. N=4 

In order to obtain the multi-particle energy spectrum, 
we have exactly diagonalized the Hamiltonian |^ nu- 
merically. As the frequency of rotation Q increases the 
ground state of the system passes through a sequence of 
states with increasing and well defined total angular mo- 
mentum L — 0,4,8,12 (see Fig. P). These states can 
be identified as follows: The state with L = is a triv- 
ial non-entangled state in which all the atoms are con- 
densed in the single particle Gaussian state with angular 
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FIG. 1: Lowest two eigenenergies (in units of hu) of the 
Hamiltonian Q for 4 particles and rj = 0.1 as a function of 
the trap rotation frequency fi/oj. The circles mark the level 
crossings and L denotes the total angular momentum of the 
ground state. The ground state sequence can be identified 
as follows (with fidelity given in brackets): L— Gaussian 
ground state (exact) , L=4 Pfaffian state (0.95), L=8 quasi- 
particle state (0.98), L=12 Laughlin state (exact). The change 
of angular momentum can readily be obtained from the in- 
creasing width of the density distribution depicted below. 



momentum m = 0. The first nontrivial ground state is 
the L = 4 state. This state is not, as one might expect, 
a single vortex state, in which all the particles would be 
condensed in the single particle state m = 1. In contrast, 
this state is highly entangled and is very close (fidelity 
0.95) to the well-known Pfaffian state: 



i<3 
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(5) 



This state is specially interesting, also in the context of 
quantum information theory, because its elementary ex- 
citations are known to exhibit non-abelian statistics . 
The next stable state in row (L = 8) can be very well 
characterized (fidelity 0.98) by a Laughlin quasiparticle 
state: 
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This state is the counterpart of the quasihole excitation, 
which has previously been studied in the context of 1/2- 
anyons in rotating Bose-Einstcin condensates . Finally, 
the last stable state is identical to the ^-Laughlin state, 
which we have already encountered in the case of two 
particles per well: 



(7) 



This state is an exact eigenstate of with zero interac- 
tion energy. In Fig. ^^e have plotted the density distri- 
bution in the xi,X2 plane of the different stable ground 



states. As the frequency of rotation rj/w increases the 
wave function spreads, and the interaction between the 
atoms decreases. 



III. ADIABATIC PATHS TO ENTANGLED 
STATES 

The sequence of entangled states we have described 
above cannot be obtained by simply adiabatically in- 
creasing the frequency of rotation il. For the rotational 
symmetry leads to level crossings between different an- 
gular momentum states (Fig. In order to pass adia- 
batically from the zero angular momentum ground state 
to higher angular momentum states the spherical sym- 
metry of the trapping potential has to be broken. For 
our optical lattice setup this can be achieved for example 
by deforming the formerly isotropic trapping potential 
on each well and letting the deformation rotate with fre- 
quency In the rotating frame the new trapping 
potential has the form Vp {uj + Aw 
the new Hamiltonian is H + H^^ with 



x{ + uj x^i and 



mO,m+2""m 



+ (m + l)aj„a,„ +h.c., (8) 



where f3m = \/(m + 2)(m + 1) and e — Auj/lu is a small 
parameter. The perturbation ||SJ leads to quadrupole 
excitations, so that states whose total angular momenta 
differ by two are coupled. 

In order to design appropriate adiabatic paths to the 
entangled states described above, we have computed nu- 
merically the energy gap between the ground and first 
excited state as a function of the parameters O/w and e 
for = 2 and iV = 4 (Fig. EJ. 

We first note that the isolines of constant energy gap 
show an approximately linear behavior. This feature can 
be easily understood from a perturbative treatment of 
the Hamiltonian ©. To first order, the energy of states 
with angular momentum L is shifted by an amount eL/4. 
Therefore the gap profile for a given e is very similar to 
the one for e = but shifted an amount '--^ e to larger 
rotation frequencies. As expected, we find that for e 7^ 
avoided crossings emerge (see Fig. The energy gap 
of the avoided crossings does, however, not in general 
increase monotonically with the deformation e. Due to 
the interplay with other excited states, "saddlepoints" 
appear in the gap profile, which makes the design of ap- 
propriate adiabatic paths a nontrivial task. For the sta- 
ble entangled states of N = 2,4 identified above these 
paths are depicted in Fig. |2 The actual time needed for 
the adiabatic path depends on the number of particles 
as well as on the state we want to achieve. For a typical 
trapping frequency to ~ (27r)30kHz and an interaction 
coupling 77 = 0.1, the evolution times for the N — 2 
Laughlin state as well as for the L = 4 and L = 8 states 
for = 4 are of the order of 10 ms. In contrast, the evo- 
lution time for the N = 4 Laughlin state is one order of 
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FIG. 2: Energy gap in units of huj between the ground and 
first excited state as a function of the rotation frequency ^/cu 
and the trap deformation e for an interaction strength rj = 0.1. 
The black lines mark appropriate paths in parameter space 
for adiabatic ground state evolution starting from the L — 
state. The adiabatic evolution times have been calculated for 
a typical trapping frequency ui ~ (27r)30kHz. Top (A'^ = 2): 
For a final fidelity = \{ip{T)\ipL)\'^ = 0.99 the Laughhn 
state (L~2) can be reached within T=6.5 ms. Bottom {N — 
4): Adiabatic path, evolution time T and fidelity T for the 
following final states (see Fig. 0: (a) Pfaffian state: T=8 ms, 
f = 0.99; (b) Quasiparticle state; T=12 ms, T = 0.99; (c) 
Laughlin state: T=215 ms, J" = 0.97. 



magnitude larger. We can understand this result in the 
following way. For the case of = 2 direct coupling of 
the L = state to the L = 2 Laughlin state is mediated 
by (jSJ. For the case of iV = 4 there is no direct coupling 
between the ground states, since their angular momenta 
differ by 4. But, as one can see from the spectrum in the 
vicinity of the crossing to the state L = 4 (Fig. ^ , there 
is a state with L = 2 near the crossing that mediates 
the coupling between the L — and the L = 4 state. 
A similar situation occurs for the crossing to the L ~ 8 
state. However, there is no such intermediate state in 
direct proximity of the crossing to the A'^ = 4 Laughlin 



state, which leads to a decrease of the energy gap by one 
order of magnitude. 

Let us also comment on the situation iV = 3. Here a 
ground state with odd angular momentum (L = 3) arises. 
From the nature of the perturbation ((SJ it is clear that 
ground state evolution is not possible. However, we have 
shown that the 1/2-Laughlin state can be reached by 
designing appropriate adiabatic paths via excited levels. 
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FIG. 3: Left side: Energy spectrum (in units hcu) for N=4 
and Tj — 0.1 in the vicinity of the first level crossing from 
the L=0 to the L=4 state (see Fig. Q left circle) . Using 
quadrupole excitations (|AL| = 2) coupling between these 
states is provided by the intermediate state L=2. Right side: 
Emergence of an avoided level crossing for a trap deformation 
e = 0.06. 



IV. FEASIBILITY 

Let us now discuss the experimental feasibility of our 
proposal for a small number of particles N. First of all 
we have assumed that the lattice wells can be treated 
independently. This requires that the overlap of the 
Wannier functions on neighbouring sites is small. We 
can estimate how intense the laser light should be in or- 
der to neglect this overlap by requiring the size of the 
highest occupied angular momentum single particle state 
(sa ^/2N — li) to be much smaller than the separation 
between lattice sites, a = ir/k. This leads to the condi- 
tion: (Vb/iJfly/'* > V^N - I/tt. Numerical calculations 
of overlap integrals between adjacent sites have confirmed 
that for N = 2(4) lattice depths of Vq/Er « 30(50) are 
required. Note that these values also guarantee the va- 
lidity of the harmonic approximation. Moreover we have 
assumed that the available single particle states on each 
well lie within the LLL. This implies that the typical 
energies per particle have to be much smaller than the 
energy gap to the next Landau level, fno. For the limit- 
ing cases of the L = state and the Laughlin state, this 
leads to the conditions (iV-l)?7/2, {N-l){l-n/uj) < 1, 
which are easily fulfilled for typical interaction strengths 
(77 ~ 0.1) and small iV. 

Finally, in order to adiabatically achieve the entan- 
gled states identified above further conditions are re- 
quired. We analyze the most restrictive case, which 
corresponds to the Laughlin state. First of all the fre- 
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quency of rotation has to be very close to the centrifu- 
gal limit. Let us find a lower bound to the critical ro- 
tation frequency at which the crossing to the Laugh- 
lin state appears. This can be done by calculating the 
rotation frequency at which the Laughlin quasiparticle 
state, tpQpiiz]) = ■ • • gf- V'L, becomes equal in en- 
ergy to the Laughlin state. §ince the quasiparticle state 
has TV units of angular momentum less than the Laugh- 
lin state and an interaction energy < 77, it follows that 
flc/uj > 1 - ri/N. For the cases of TV = 2(4) this condi- 
tion is in agreement with the exact values found above. 
Secondly, the evolution time required for the adiabatic 
path has to be much smaller than the typical decoher- 
ence time. We can estimate this time in the following 
way. Given the critical frequency above and that the 
position of the avoided crossing is displaced to larger ro- 
tation frequencies an amount proportional to e it follows 
that the maximum e wc can have is ~ t]/N, corresponding 
to a rotation frequency Q/uj — 1. Assuming an energy 
gap « e it follows that the typical evolution time scales as 
Nr]. For the case of A'' = 2(4) and typical rj and uj these 
times are of the order of tens of miliseconds as exactly 
found above, which is much smaller than the typical life 
time of the lattice states. Finally, a high degree of con- 
trol of the parameters il/uj and e is required to perform 
the appropriate adiabatic paths. The required precision 
scales again as rj/N, which for the case of A = 4 means 
a control of the parameter space up to the second digit. 

From our analysis it follows that the adiabatic creation 
of the Laughlin state by means of low angular momentum 
excitations, as quadrupole excitations, becomes very dif- 
ficult in samples with large number of particles jlTl [l8| . 
Even if the centrifugal limit is possible to achieve, as it 
happens when including an additional trapping po- 
tential the adiabatic creation of the Laughlin state 
is still very demanding. One reason is that the rotation 
frequency and the trap deformation have to be controlled 
within a precision that also scales linearly with A^. Fur- 
thermore we point out that only the exact knowledge of 
the multi-particle energy spectrum allows to design adi- 
abatic paths that minimize the evolution time. 



V. DETECTION 

We discuss now how the entangled states described 
above can be detected experimentally by measuring dif- 
ferent properties of the states. As an important feature 
of our lattice setup of independent wells, we note that 
any signal will be highly enhanced by a factor equal to 
the number of occupied lattice sites ( ^ 150, 000 Q). 

i) Density profiles. A very characteristic feature of the 
entangled states with large angular momentum that we 
have described is that they exhibit a much more extended 
density distribution than the non-entangled L = state, 
in which the particles are much more confined in space. 
For the 1/2-Laughlin state the typical radius is given by 
f « \/2A - 1 £. In the case of A = 2(4) this results in 



a radius that is ~ 2(3) times lar ger than in the case of 
the condensate. As proposed in [l9j the density profile 
of states within the LLL can be measured in a time of 
flight (TOF) image of the atomic system, since the mo- 
mentum distribution coincides with the density profile for 
LLL states. In our case of independent 3D wells, a TOF 
absorption picture after expansion time t will exhibit a 
broad central peak of the form: 

^ (^ \poi-iz/iu;t),X3/icjt)\^ . (9) 

Here, po{z, X3) is the initial density distribution on a sin- 
gle well. In the TOF image it is enhanced by a factor pro- 
portional to the number of lattice sites A^^ and rescaled 
by a factor ujt ^ 1. The tt/2 rotation z —iz leaves 
isotropic states, like the FQH states described above, un- 
affected. The underlying assumption of free (interaction- 
less) expansion is justified, since the interaction energy is 
small compared to the kinetic energy (in the stationary 
frame) . 

i\)Angular momentum. For any state within the 
LLL integration over the density distribution gives 
J dr r^p(r) = L + N . Thus in the limit of weak in- 
teraction the total angular momentum can be extracted 
directly from the TOF picture. 

iii) Correlation functions. Here we propose a novel 
technique that makes directly use of the rich possibil- 
ities offered by the optical lattice setup and which al- 
lows to measure both the gi = {iP\y)iIj{y')) and 52 — 
{ip''' {r)ip''' {r')ip{r)ip{r')) correlation functions. The 52 cor- 
relation function is for instance very characteristic for 
a Laughlin state. Since particles can only be at least 
in relative angular momentum = 2 it follows that 
.92 c>c |r — r'\'^. This behavior reveals the ^-fractional 
nature of this Laughlin state. 

We consider two species a and b (hyperfine levels) of 
bosonic atoms, which can be coupled via Raman transi- 
tions. We start with atoms in level a and create the en- 
tangled state of interest with the method described 
above. Next we apply a 7r/2-pulse with the laser and 
create an equal superposition of a and b states. Finally, 
we shift the lattice potential trapping atoms of type b 
(as proposed in and realized in ^) by a distance 
To small compared to the lattice spacing and perform 
another 7r/2-pulse. In the Heisenberg picture this proce- 
dure corresponds to the following transformation of the 
field operator for species a: 

i^a(r)-^ V,(r)+V,(r + ro) . (10) 

Thus the density distribution of atoms of type a in this 
new state contains information about the gi corre- 
lation function of the original state: 

{^Mir)Mr)\^f) - (11) 
(^-.K^t (r) + (r + ro))(V^a(r) + + ro))|*.) • 

Using this procedure we can also measure higher order 
correlation functions like g2. In this case measuring the 
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interaction energy of the final state will allow us to cal- 
culate the g2 of the initial state. For instance, for the 
Laughlin state we have: 

E^nt{ro) = (12) 
^ j dr(*,|^l(r)Vl(r + ro)V'a(r)V'a(r + ro)|*,) . 

The interaction energy is, unfortunately, not directly ac- 
cessible experimentally. However, the total energy of the 
final state can be obtained from integrating over the TOF 
absorption picture, since energy is conserved during the 
time of flight. For small coupling -q, however, the mea- 
surable effect due to interactions will be small compared 
to the kinetic part of the energy. In addition, the ki- 
netic energy itself shows a significant dependence on the 
shifting ro, which has to be distinguished from the inter- 
action. Hence, we propose to tune the scattering length 
as (e.g. via a photo association induced Feshbach reso- 
nance |2lj |) and to measure the interaction energy both 
in the weak and strong scattering limit. The difference 
would then reveal the characteristic behavior of the 32 
correlation function. 



We finally note that, as a further way of detection for 
the = 4 Laughlin state, a strong reduction of the three 
body losses should be observed. 



VI. CONCLUSION 

In conclusion, we have shown how to motionally entan- 
gle a small number of particles into a sequence of inter- 
esting FQH states. We have fully addressed the adiabatic 
creation of these states and proposed new techniques for 
their experimental detection. 
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